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Maucha 图 的 Matlab 实现 及 其 在 水 化 学 中 的 应 用 
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摘 要 : 水 中 可 溶性 离子 的 含量 及 各 离子 的 比例 关系 是 衡量 水 质 的 重要 指标 ,但 常用 的 方法 不 能 同 
时 反映 这 两 种 关系 。Maucha 图 既 可 以 反映 8 种 主要 阴 、 阳 离子 的 总 浓度 ,同时 也 可 以 指示 各 离子 
的 相对 浓度 ,兼顾 了 离子 的 绝对 含量 和 相对 含量 ,具有 较 强 的 实用 性 。 但 由 于 缺乏 绘制 Maucha 图 


的 专业 软件 ,导致 Maucha 图 在 实际 应 用 中 受到 


的 分 析 提 供 和 参考 依据 。 
关 键 词 : Maucha 图 ; 水 化 学 ; Matlab; 绘图 ; 应 用 


水 资源 是 人 类 赖 以 生存 和 发 展 的 基础 ,但 我 国 
水 资源 时 空 分 布 极 不 均匀 、 人 均 水 资源 占有 量 少 , 水 
资源 安全 面临 很 大 挑战 ,特别 是 部 分 天 然 水 体 
水 质 不 达标 ,无 法 满足 人 们 生产 和 生活 的 需要 ,因而 
水 质 分 析 成 为 水 资源 开发 利用 的 前 提 之 一 。 各 类 水 
化 学 图 被 广泛 应 用 于 水 质 分 析 ” ,可 反映 水 中 可 
溶性 离子 的 含量 及 其 比例 关系 ,其 中 常用 折线 图 来 
反映 水 中 各 离子 的 绝对 含量 。 该 图 以 横 坐 标 表示 各 
类 离子 , 纵 坐 标 表 示 离 子 浓度 。 折 线 图 虽然 能 准确 
反映 各 离子 的 真实 浓度 ,但 当 同 一 个 折线 图 内 样品 
量 较 多 时 ,数据 点 之 间 的 重合 可 能 导致 无 法 有 效 区 
分 样 点 。 此 外 ,折线 图 反映 的 是 离子 的 绝对 浓度 ,并 
不 能 反映 单个 离子 占 所 有 阴 / 阳 离子 或 总 离子 的 比 
例 , 即 不 能 反映 单个 离子 的 相对 浓度 。 

Maucha 图 很 好 的 解决 了 上 述 问 题 , 既 可 以 避免 
样品 间 的 互相 重 铸 ,又 可 以 指示 样品 各 离子 的 相对 
浓度 ,其 大 小 还 可 以 反映 离子 的 绝对 浓度 和 矿 化 度 。 
该 图 由 匈牙利 科学 家 Maucha 于 1932 年 以 德语 著作 
的 形式 提出 " ,后 经 其 他 科学 家 数 次 修改 。 例 如 ， 
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很 大 限制 。 通 过 介绍 Maucha 图 的 演化 过 程 , 绘 条 
Maucha 图 时 建立 坐标 系 及 确定 Maucha 图 各 组 成 部 分 坐标 值 的 方法 ,并 阅 述 了 基于 Matlab 语言 实 
Il Maucha 图 计算 机 输出 的 编程 流程 。 最 后 以 内 蒙古 浑 善 达 克 沙 地 东部 地 区 为 例 , 说 明了 Maucha 
图 在 天 然 水 体 水 化 学 分 析 中 的 应 用 。 该 应 用 表明 ,Maucha 图 在 反映 样品 矿 化 度 、 指 示 各 离子 的 相 
对 浓度 和 绝对 浓度 、 判 断水 体 水 化 学 类 型 及 水 化 学 数据 空间 分 析 等 方面 均 具 有 和 较 强 的 优势 ,能 显著 
提高 水 化 学 数据 的 可 视 化 表达 能 力 ,丰富 水 化 学 图 
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在 水 文 地 质 中 的 应 用 ,并 为 水 化 学 数据 更 深层 次 


原始 型 Maucha 图 只 能 指示 各 离子 的 相对 浓度 ,不 
能 反映 样品 的 矿 化 度 ,因此 HEDCPETH 在 原始 型 
Maucha 图 旁 增加 了 一 个 对 数 标尺 ,用 以 反映 样品 矿 
化 度 的 大 小 ,并 对 矿 化 度 小 于 等 于 海水 的 样品 用 虚 
线 绘制 ,而 对 矿 化 度 大 于 海水 的 样品 用 实 线 绘制 。 
但 该 改进 方法 仍 无 法 直观 的 反映 不 同样 品 矿 化 度 的 
差异 。 基 于 这 个 问题 ,BROCH'”' 为 Maucha 图 的 圆 
形 部 分 引入 了 内 切 正 十 六 边 形 的 概念 ,该 内 切 正 十 
六 边 形 并 不 需要 绘制 在 Maucha 图 中 ,但 可 以 用 其 
面积 来 表示 样品 中 8 个 主要 阴 、 阳 离子 的 总 浓度 。 
由 于 这 8 大 离子 是 产生 水 样 矿 化 度 的 主要 成 分 ， 
而 可 以 用 其 近似 反映 样品 的 矿 化 度 。 改 进 后 的 
Maucha 图 既 可 以 通过 内 切 正 十 六 边 形 的 面积 来 反 
映 样 品 的 绝对 含量 ,又 可 以 通过 表示 各 离子 含量 的 
四 边 形 的 面积 指示 其 相对 含量 。 因 圆 的 面积 与 其 内 
切 正 十 六 边 形 的 面积 相对 应 ,因此 实际 应 用 中 可 直 
接 根据 圆 的 面积 来 比较 不 同样 品 的 矿 化 度 。 本 文 以 
下 所 述 均 指 经 BROCH 改进 后 的 Maucha 图 。 
Maucha 图 与 雷达 图 (Radarchart ) 在 外 观 上 较 
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任 孝 宗 等 :Maucha 图 的 Matlab 实现 及 其 在 水 化 学 中 的 应 用 


为 相似 ,但 Maucha 图 并 不 能 像 雷 达 图 那样 用 常见 
科学 绘图 软件 实现 。HASSELL 基于 Qbasic 语言 实 
IT Maucha 图 的 计算 机 输出 中 ,但 基于 Qbasic 的 
操作 系统 过 于 陈旧 ,导致 该 方法 不 能 得 到 有 效 推广 。 
Matlab 语言 可 在 当前 多 数 主流 计算 机 操作 系统 上 运 
行 , 是 近年 来 被 广大 科技 工作 者 广泛 使 用 的 计算 机 
编程 语言 ,其 不 但 具有 可 读 性 强 .调试 简单 等 特点 ， 
且 其 强大 的 绘图 功能 在 水 化 学 图 的 计算 机 输出 方面 
表现 出 很 大 优势 "1 。 本 文 在 阐明 Maucha 图 绘制 原 
理 的 基础 上 ,以 Matlab 语言 编程 实现 了 Maucha 网 
的 计算 机 输出 ,并 通过 实例 说 明 Maucha 图 在 水 化 
学 中 的 应 用 ,以 期 扩展 该 图 在 水 化 学 中 的 应 用 范围 。 


1 Maucha 图 的 组 成 及 其 含义 
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代表 其 他 离子 的 另外 7 个 四 边 形 , 且 这 8 个 四 边 形 
的 面积 之 和 等 于 圆 内 切 正 十 六 边 形 的 面积 。 需 要 注 
意 的 是 ,Maucha 图 中 的 离子 浓度 指 的 是 当量 浓度 ， 
且 阴 、 阳 离子 电荷 数 相等 , 呈 电 中 性 。 


2 Maucha 图 的 绘制 原理 


(1) 单位 转换 。 离 子 浓度 一 般 以 mg + LR 
示 , 因 此 需要 将 其 转换 为 当量 浓度 (meq ' 工 ” ) 。 

(2) 阴 、 阳 离子 电荷 平衡 及 质量 平衡 原理 检验 。 
Maucha 图 用 两 个 面积 相等 的 半圆 分 别 指示 阴 、 阳 离 
子 总 浓度 ,其 前 提 条 件 是 阴 、 阳 离子 总 电荷 数 需 符 合 
电荷 平衡 原理 。 这 是 衡量 水 化 学 数据 准确 性 和 精度 
的 标准 ,同时 也 是 能 否 应 用 Maucha 图 的 先决 条 件 。 


= 因此 在 绘制 Maucha 图 前 需要 对 阴 、 阳 离子 电荷 平 
in Maucha 图 由 1 个 加 和 主体 位 于 圆 内 的 8 个 四 衡 进 行 检验 。 由 于 水 中 除 Maucha 图 指示 的 八 种 主 


LO 边 形 组 成 , 圆 又 被 分 为 左右 两 个 半圆 。 左 半圆 由 四 ”要 离子 外 还 含有 其 他 带电 微量 离子 .胶体 以 及 有 机 


A 个 四 边 形 组 成 ,其 面积 分 别 指示 四 种 阴离子 ( 按 逆 


时 针 顺 序 从 上 到 下 依次 为 C03-、HCO; 、CL 和 
S04 ) 的 浓度 ;相对 应 右 半 圆 也 由 四 个 四 边 形 组 成 ， 
其 面积 分 别 指示 四 种 阳离子 ( 按 逆 时 针 顺 序 从 下 到 
上 依次 为 Mg *、Ca*、Na’' AK" ) 的 浓度 (图 1)。 
圆 的 半径 大 小 由 其 内 切 正 十 六 边 形 的 面积 决定 ,其 
内 切 正 十 六 边 形 的 面积 如 前 文 所 述 由 8 大 离子 的 总 

代表 某 个 离子 的 四 边 形 的 确定 方法 为 : 当 圆 被 
等 分 为 8 个 悄 形 后 ,对 于 每 个 肩 形 而 言 ,由 组 成 鹿 形 
的 两 条 半径 及 扇形 内 或 外 的 一 个 点 (该 点 须 位 于 过 
圆 的 十 六 等 分 点 的 射线 上 ) 可 确定 一 个 四 边 形 ,该 
四 边 形 的 面积 须 等 于 该 离子 的 浓度 。 同 理 , 可 确定 


阳离子 
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| 

l 


图 1 Maucha 图 及 其 含义 


Fig.1 Maucha diagram and its implications 


物 等 ,因此 阴 、 阳 离子 电荷 平衡 常数 满足 + 上 5% 时 通 
常 被 认为 是 符合 电荷 平衡 原理 的 。 电 荷 平衡 党 
BEE. B. 的 计算 公式 "2 如 下 所 示 : 
_SC+SA 
SC - SA 
式 中 :SC 为 阳离子 电荷 总 数 ; SA 为 阴离子 电荷 总 
数 。 此 外 ,水 化 学 数据 还 需 满足 质量 平衡 原理 ,要 求 
样品 矿 化 度 等 于 8 大 离子 含量 之 和 , 旦 精度 控制 在 
+5% 。 

(3) 建立 坐标 系 ,确定 坐标 原点 。 如 图 2 和 图 
3, 以 圆心 0 为 坐标 原点 ,0G 方向 为 x Hh, OD 方向 
H y 轴 , 建 立 坐 标 系 。 

(4) 圆 的 半径 R 和 顶点 B 的 确定 。 当 COS  、 
HCO; .Cl .S01 Mg”*、Ca' Na 和’ 中 某 离子 
的 含量 小 于 8 大 离子 平均 值 时 ,代表 某 离 子 含量 的 
四 边 形 将 位 于 圆 内 (图 2) ;相应 的 , 当 某 离子 的 含量 
大 于 8 大 离子 平均 值 时 ,代表 某 离 子 含量 的 四 边 形 
将 位 于 圆 外 (图 3 ) ; 当 某 离子 的 含量 等 于 8 大 离子 
平均 值 时 ,代表 某 离 子 含量 的 四 边 形 将 位 于 圆 上 。 
无 论 上 述 哪 种 情况 ,在 确定 Maucha 图 各 参数 的 过 
程 中 均 使 用 相同 的 公式 ,因此 ,这 里 以 代表 KK’ 的 四 
边 形 ODBC 为 例 说 明 圆 的 半径 R 和 四 边 形 顶 点 B 
的 确定 过 程 。 图 2 和 图 3 中 D.C.G.HIJK 和 L 
为 圆 的 8 个 八 等 分 点 ,A 为 圆 的 1 个 十 六 等 分 点 , 线 
Be AC 为 圆 内 切 正 十 六 边 形 的 一 条 边 ,线段 CF 垂直 


E. B. 


x 100% (1) 
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K2 Maucha 图 各 组 成 部 分 坐标 的 确定 方法 
(顶点 B 位 于 圆 内 ) 


Fig.2 Determining the coordinates of the components in 


Maucha diagram ( vertex B inside the circle) 


图 3 Maucha 图 各 组 成 部 分 坐标 的 确定 方法 
(顶点 B 位 于 圆 外 ) 


Fig.3 Determining the coordinates of the components in 


Maucha diagram ( vertex B outside the circle ) 


于 线段 OA, 线段 AE 垂直 于 线段 OC, ZAOC 为 
22.5。。 
根据 Maucha 图 的 组 成 特点 可 知 : 圆 内 切 正 十 
六 边 形 的 面积 等 于 16 倍 的 三 角形 OAC 的 面积 ,也 
等 于 8 大 离子 当量 浓度 的 总 和 ,结合 三 角 函 数 关系 
可 得 : 
R? . sin 22.5°/2 = T/16 (2) 


式 中 :及 为 圆 的 半径 ;了 为 8 大 离子 的 当量 浓度 总 
和 。 
三 角形 OBC 的 面积 等 于 K* 当量 浓度 的 一 半 ， 
结合 三 角 函 数 关系 可 得 : 
b + R - sin 22.5°/2 =P/2 (3) 


式 中 :2 为 线段 OB 的 长 度 ;P 为 K* 当量 浓度 。 
顶点 B 的 坐标 为 : 


bx =b + sin 22.5° (4) 
by =b + cos 22.5° (5) 
式 中 :bx by 分 别 为 顶点 B 的 横 、 纵 坐标 。 


由 公式 (3) 求 得 ,再 带 入 (4) 和 (5 ) 即 可 求 得 
顶点 B 的 坐标 bx by. 

同 理 ,另外 7 个 四 边 形 的 顶点 也 可 以 由 上 述 方 
法 确定 。 


3 Maucha 图 的 编程 流程 


Maucha 图 的 Matlab 实现 流程 见 图 4。 首 先 , 将 
被 测 样品 的 离子 浓度 (mg L” ) 转化 为 当量 浓度 
(meq . 工 - )。 其 次 ,判断 样品 是 否 符合 电荷 平衡 原 
理 及 质量 平衡 原理 。 如 果 电 蓓 不 平衡 或 不 满足 质量 
平衡 原理 , 则 在 程序 运行 窗口 显示 :“ 样 品 电 荷 不 平 
衡 或 不 满足 质量 平衡 原理 ,无 法 绘制 Maucha RI”, 
然后 退出 程序 ;如 果 阴 、 阳 离子 电荷 数 平 衔 且 满足 质 
量 平 衔 原理 ,程序 继续 运行 。 再 次 ,以 圆心 为 坐标 原 
点 \ 两 条 相互 垂直 的 半径 方向 为 x 轴 和 y 轴 建 立 直 
角 坐 标 系 。 根 据 上 文 所 述 公 式 (1) 计 算 圆 的 半径 
及 ,并 用 plot 命令 绘制 圆 ,然后 根据 公式 (3) (4) 和 
(5) 及 循环 语句 计算 顶点 了 B 的 坐标 bx by 以 及 另外 
七 个 四 边 形 的 顶点 。 最 后 用 plot 命令 绘制 出 8 个 四 


边 形 。 
数据 输入 (mg + L) 


转 | 换 
HEHE meq- L) 


确定 坐标 原点 并 
建立 坐标 系统 


判断 样品 是 否 符合 电荷 
平衡 原理 和 质量 平衡 原理 
v 
计算 半径 R 及 顶点 B 
的 坐标 bx、by 
> < 
绘制 输出 Maucha 图 


图 4 Maucha 图 的 编程 流程 
Fig.4 Program flow for plotting Maucha Diagram 
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在 整个 Maucha 图 的 Matlab 实现 过 程 中 ,几乎 
不 出 现 精度 的 损失 。 首 先 ,将 离子 浓度 从 mg * LL” 
转换 为 当量 浓度 meq ' 工 ”时 ,所 采用 的 原子 量 数值 
取 3 位 有 效 数字 ,单位 转换 过 程 中 的 数值 类 型 采用 
Matlab 中 精度 最 高 的 双 精 度 浮 点 数 。 其 次 ,在 计算 
圆 的 半径 只 和 顶点 B 的 坐标 时 ,数值 类 型 仍然 采用 
双 精 度 浮 点 数 。 以 上 措施 能 最 大 限度 的 保障 结果 具 
有 和 较 高 精度 。 


4 ” Maucha 图 在 浑 善 达 克 沙 地 的 应 用 


Maucha 图 在 水 化 学 分 析 中 使 用 较为 广泛 ,其 应 
用 方式 可 归纳 为 两 种 :一 是 将 不 同样 品 的 Maucha 
图 绘制 在 同一 图 中 ,以 便 进行 样品 间 的 比较 ,例如 
HAROLD 等 将 Maucha 图 应 用 于 南非 的 研究 工 
Ye?! ;其 二 是 将 样品 的 Maucha 图 绘制 到 研究 区 概 
况 图 中 ,利用 Maucha 图 对 水 化 学 数据 进行 空间 分 
析 , 例 如 CONZONNO 和 ULIBARRENA 将 Maucha 图 
应 用 于 阿根廷 的 研究 工作 "1。 

本 文 以 浑 善 达 克 沙 地 天 然 水 体 为 例 ,采用 公开 
发 表 的 水 化 学 数据 ,说明 Maucha 图 在 水 化 学 中 
的 应 用 。 浑 善 达 克 沙 地 位 于 我 国内 蒙古 高 原 东部 ， 
主要 以 固定 和 半 固 定 沙丘 为 主 “”i。 沙 地 天 然 水 
体 的 COS” 含量 很 低 或 者 不 含 C03 ,而 其 天 然 水 体 
特别 是 地 下 水 由 于 受到 人 类 活动 的 影响 NOS 含量 
较 高 ,研究 NO, 含量 的 变化 比 研究 CO; -含量 的 变 
化 更 有 意义 。 因 此 本 文 对 Maucha 图 做 了 调整 ,用 
NO; 替代 CO; ,并 在 等 比例 条 件 下 将 浑 善 达 克 沙 
地 东部 地 区 天 然 水 体 的 Maucha 图 绘制 在 同一 图 中 
(图 5)。 绘 制 Maucha 图 前 , 先 对 其 是 否 满足 电荷 平 
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图 5 Maucha 图 在 浑 善 达 克 沙 地 的 应 用 
Fig.5 Application of Maucha diagram in Hunshandake 
Sandy Land 
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衡 原理 及 质量 平衡 原理 进行 判断 ,经 程序 判断 ,该 数 
据 满足 电荷 平衡 原理 及 质量 平衡 原理 。 图 中 gl 、g2 
aaa g11 为 地 下 水 样品 , 共 11 个 ; 1 .2……]16 为 湖 


sl 和 s2 为 果 水 样品 , 共 2 个 。 

如 前 文 所 述 ,可 用 圆 的 面积 指示 样品 的 矿 化 度 。 
从 图 5 可 知 ,大 多 数 地 下 水 样品 的 矿 化 度 高 于 湖水 、 
河水 和 泉水 样品 ,这 可 能 和 地 下 水 与 含水 层 物质 的 
水 一 岩 作用 时 间 长 溶解 较 多 的 物质 有 关 …] 。 和 典 
型 西部 沙漠 如 巴 丹 吉林 沙漠 的 湖水 样品 相 比 ， 
浑 善 达 克 沙 地 东部 地 区 湖水 样品 的 矿 化 度 并 不 高 ， 
暗示 该 地 区 的 共 发 量 小 于 我 国 西部 沙漠 。 通 过 对 比 
代表 各 离子 浓度 的 四 边 形 的 面积 可 以 判断 该 样品 的 
优势 阴 、 阳 离子 。 如 图 5 所 示 ,地 下 水 各 样品 的 优势 
阴离子 并 不 相同 ,为 NO; 或 HCO; BR S0} 或 Cl ; 
而 大 多 数 地 下 水 样品 的 优势 阳离子 较为 一 致 , 主要 
是 Cat ,其 次 是 Na' 。 湖 水 、 河 水 和 泉水 样品 的 优 
势 阴 .阳离子 相对 一 致 , 除 凡 水 样品 s1 (优势 阴离子 
为 S04 ,其 次 为 HCO; ) 外 ,其 他 样品 的 优势 阴离子 
HJ HCO, ,而 优势 阳离子 和 地 下 水 的 优势 阳离子 
相同 ,主要 是 Ca ,其 次 是 Na *。 

优势 离子 也 可 以 通过 计算 获取 ,一 般 将 百 分 含 
量 (基于 摩尔 浓度 ) 大 于 25% 的 阴 、 阳 离子 统计 制 成 
水 化 学 类 型 表 '“1 。 结 果 表 明 ( 表 1) ,多 数 地 下 水 样 
品 的 优势 阴离子 较 复 杂 ,而 优势 阳离子 主要 是 Ca 
和 Na* ;多数 地 表 水 的 优势 阴离子 是 HCO， 而 优势 
阳离子 是 Ca 和 Na* 。 该 结果 与 基于 Maucha 图 的 
分 析 基 本 一 致 但 略 有 差别 (原因 是 Maucha 图 使 用 
当量 浓度 而 水 化 学 类 型 表 使 用 摩尔 浓度 ) ,表明 
Maucha 图 可 以 直观 的 反映 样品 的 优势 离子 。 

将 不 同样 品 的 Maucha 图 按 等 比例 原则 绘制 在 
采样 图 的 对 应 位 置 ,可 用 以 分 析 样 品 的 矿 化 度 及 各 
离子 相对 浓度 在 空间 上 的 变化 关系 。 本 文选 取 了 9 
个 样品 ,将 其 绘制 在 DEM 底 图 上 (图 6)。 从 图 6 中 
可 以 看 到 ,地 下 水 样品 主要 分 布 在 西 拉 木 伦 河 源头 
西南 方向 ,其 矿 化 度 高 于 地 表 水 样品 ,其 与 其 他 地 表 
水 样品 的 距离 约 为 50 km ,因此 气候 并 不 是 引起 它 
们 较 大 差异 的 主要 因素 。 通 过 氧 氧 同位 素 等 指标 的 
进一步 分 析 ,发 现 引起 这 种 差异 的 主要 因素 是 水 体 
的 补给 方式 不 同 中 | 。 

该 空间 分 析 方 法 存在 一 定局 限 性 ,例如 ,由 于 
Maucha 图 所 占 面积 较 大 , 当 样 品 较 多 时 ,很 难 将 所 
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Rl 浑 善 达 克 沙 地 天 然 水 体 水 化 学 分 类 表 
Tab.1 Hydrochemical classification of natural water in 
Hunshandake Sandy Land 


水 体 


样品 
编号 类 型 水 化 学 类 型 
gl 也 下 水 HCO; -Ca 
g2 也 下 水  HCOș - Ca2+ 
3 也 下 水 HCO;y -C17 -Ca?* -Nat 
g4 也 下 水 Cl- -HCO; -S02?- -Ca?* -Na* -Mg + 
25 也 下 水 Cl- -HCO; -S03?- -Ca2+ -Na+ 
26 也 下 水 Cl- -HCO; - Na* - Ca2+ 
27 也 下 水 NO, 一 Ca2+ - Nat 
28 WFK HCO; -NO; -Na+ -Ca2+ 
29 岂 下 水 HC03 -Ca?* 
g10 也 下 水 HCO; - NO; -Ca?* -Na+ 
ell ih F7k HCO; - SO -Ca?* - Nat 
ll 湖水 HCO; - Na* - Ca2+ 
湖水 HCO; - Na* - Ca2+ 
B 湖水 HCO; -Ca —Na* - Mg2+ 
4 湖水 HCO; - K* - Ca2+ - Mg?* 
15 湖水 HCO; - Na* - Ca2+ 
16 湖水 HCO; - Na* - Ca** - Mgt 
rl 河水 HCO; - Na* - Ca2+ 
12 河水 HCO; - Na* - Ca2+ 
3 河水 HCO; - Na* - Ca2+ 
4 河水 HCO; - Ca?* - Na* 
15 湖水 HCO; - Ca?* - Na+ 
sl 泉水 HCO; - S07- - Na* - Ca?* 
s2 泉水 HCO; - Ca?* - Na+ 
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图 6 ”Maucha 图 在 水 化 学 数据 空间 分 析 中 的 应 用 


Fig.6 Application of Maucha diagram in spatial analysis 


for hydrochemical data 


有 样品 的 Maucha 图 绘制 到 一 幅 底 图 上 ,限制 了 
Maucha 图 在 空间 分 析 中 的 应 用 。 针 对 应 用 出 现 的 
问题 ,有 两 种 解决 办 法 。 其 一 ,绘制 大 比例 尺 的 研究 
区 底 图 ,其 二 , 按 比 例 缩小 Maucha 图 ,以 便 最 大 限 
度 的 将 所 有 样品 的 Maucha 图 涵盖 到 同一 研究 区 底 
图 上 。 另 外 ,也 可 以 通过 将 所 有 样品 的 Maucha 图 
绘制 到 一 幅 图 中 ,然后 对 照 研 究 区 概况 图 ,间接 实现 
Maucha 图 的 空间 分 析 功 能 。 


5 结论 . 


Maucha 图 既 可 以 反映 8 大 离子 的 矿 化 度 ,又 可 
以 指示 各 离子 的 相对 浓度 和 绝对 浓度 ,具有 较 强 的 
实用 性 ,但 由 于 缺乏 绘制 Maucha 图 的 专业 软件 导 
致 该 方法 不 能 得 到 普遍 应 用 。 本 文 介绍 了 Maucha 
图 的 演化 过 程 ,阐述 了 绘制 Maucha 图 时 坐标 系统 
的 建立 及 其 各 组 成 部 分 ( 圆 半径 和 四 边 形 顶点 ) 坐 
标 值 的 确定 方法 ,并 以 Matlab 语言 为 工具 说 明了 实 
现 计算 机 输出 Maucha 图 的 编程 流程 。 在 此 基础 
上 ,本 文 以 内 蒙古 浑 善 达 克 沙 地 东部 地 区 天 然 水 体 
水 化 学 数据 为 例 ,根据 研究 区 天 然 水 体 水 化 学 特点 ， 
用 NO; 替换 C0; ”绘制 了 基于 Matlab 语言 的 
Maucha 图 。 结 果 表 明 , Maucha 图 在 反映 样品 矿 化 
E .指示 各 离子 相对 浓度 和 绝对 浓度 .判断 水 体 水 化 
学 类 型 及 水 化 学 数据 空间 分 析 等 方面 具有 较 强 的 优 
势 。 相 比 其 他 编程 语言 ,基于 Matlab 绘制 的 Maucha 
图 具有 更 广泛 的 计算 机 操作 系统 基础 .更 好 的 可 读 
性 和 更 强 的 可 视 化 能 力 。 此 外 ,在 该 方法 的 基础 上 
还 可 以 进行 二 次 开发 ,根据 研究 区 的 水 化 学 特点 ,用 
其 他 离子 蔡 换 原 Maucha 图 中 的 离子 (如 用 NO; 替 
换 CO; ) ,以 实现 特定 的 研究 目的 。 该 方法 的 应 用 
将 有 助 于 推动 Maucha 图 在 水 文 地 球 化 学 过 程 的 识 
别 和 解释 中 的 应 用 ,提高 水 化 学 数据 的 可 视 化 能 
为 更 深层 次 的 水 化 学 数据 分 析 提 供 参 考 依据 。 
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Plotting Maucha Diagram based on Matlab and its applications 
on hydrochemistry 


REN Xiao-zong', LIU Min', LI Jian-gang , LI Ji-yan' 
(1 Taiyuan Normal University ,School of Geographical Sciences , Jinzhong 030619 , Shanxi , China ; 


2 The construction and management station of three North Shelterbelt Project in Baiyin , Baiyin 730900 , Gansu , China) 


Abstract: The content and the proportion of soluble ions in water are important indicators to assess water quality , 
but common methods used in hydrochemistry cannot reflect the two relationships simultaneously. Maucha diagram 


,HCO; ‚C17 ,SO ， 


not only reflects the total concentration of eight main anions and cations in water (i. e. CO3; 


Mg** ,Ca ,Na* and K* ) ,but also indicates relative and absolute concentrations for each ion simultaneously. 
Thus ,the Maucha diagram is deemed to be very practical. However, due to the lack of professional software to plot 
Maucha diagram, the application of this diagram is very limited in hydrochemical research. In this paper, firstly , we 
introduced the evolution history of Maucha diagram which was proposed by Maucha in 1932 and developed by 
Hedgpeth and BROCH. Secondly ,we introduced the components of Maucha diagram and its implications. This dia- 
gram makes up of a circle and eight quadrangles inside. The area of the tangent 16-sided polygon in the circle is uti- 
lized to represent the total concentration of eight anions and cations, and the area of the eight quadrangles represents 
the content of eight ions. In addition ,the total cations and anions are supposed to submit the charge balance princi- 
ple and the mass balance principle. Thirdly ,the methods of how to establish a coordinate system and determine the 
coordinates of the circle and quadrangles in Maucha diagram were introduced in detail. Fourthly , the program flow of 
how to plot Maucha diagram based on Matlab was stated. Finally , taking the eastern Hunshandake Sandy Land in In- 
ner Mongolia, China as an example ,this paper illustrated the applications of Maucha diagram in the field of hydro- 
chemical analysis for natural waters. The Hunshandake Sandy Land belongs to a climate transition zone near the 
middle-latitude boundary of semi-humid and semi-arid to arid climate. Elevations in Hunshandake Sandy Land de- 
crease from ca. 1 300 m in the southeast to ca. 1 000 m a. s. l. in the northwest, while mean annual precipitation de- 
creasing from ~450 mm in the southeast to 150 mm or so in the northwest. Since natural water in Hunshandake 
Sandy Land is characteristic of high concentration of NO; and very low concentration of CO3” ,we used NO} in- 
stead of CO; to plot Maucha diagram with Matlab. Before plotted Maucha diagram ,we checked the ion data and on- 
ly charge balance and mass balance data were used. Applications show that Maucha diagram has strong advantages 
in reflecting the total concentration of anions and cations, indicating the relative and absolute concentrations ions, 
judging the hydrochemical type of water body and spatial analysis for hydrochemical data. Compared to other com- 
puter languages ,the Maucha diagram program written by Matlab can run on almost all kinds of computer systems , 
and have better readability and visualization. Moreover, based on hydrochemical characteristics of certain study are- 
a,the program can be developed to adapt to new situations,such as in this paper we use NO; rather than CO3 to 
plot Maucha diagram. Thus , the application of the diagram may improve the visual expression ability for hydrochemi- 
cal data, enrich the application of diagrams in hydrochemistry , and provide a reference for further analysis of hydro- 
chemical data. 


Key words: Maucha Diagram; hydrochemistry; Matlab; plotting; applications 


